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摘 要 : 由 于 考虑 系统 的 安全 性 和 经 济 性 无 法 进行 开 环 实验 ， 而 目前 最 常用 的 闭环 继 电 反 馈 辩 识 方法 需要 进行 数 次 的 
继 电 反馈 实验 ， 并 且 该 方法 只 能 状 识 国定 结构 的 低 阶 模型 ， 无 法 满足 实际 工业 的 需要 。 针 对 这 一 情况 ， 提 出 了 一 种 新 
的 含 积分 环节 过 程 对 象 的 频 域 模型 辨识 方法 ， 其 能 快速 找到 对 象 模型 的 主要 频率 段 ， 分 析 频 率 段 内 的 频率 特性 ， 从 而 
准确 的 辨识 出 含 积 分 环节 对 象 的 模型 。 仿 真实 验 的 结果 表明 该 辨识 方法 能 够 高 效 且 准 确 的 辨识 出 含 积 分 加 环节 的 对 象 。 
关键 词 : 模型 辨识 ; 含 积分 环节 ; 频率 响应 

中 图 分 类 号 : TP391.4 doi: 10.3969/j.issn.1001-3695.2018.02.0089 


Modeling identification for process with integral unit in frequency domain 
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Abstract: The model with integral unit is one process object in actual industry especially in chemical engineering and thermal 
power plant. The open loop experiment could not be done because of systems’” security and economy. Also the close-loop 
identification method which is mostly commonly used to identify the object is by using relay feedback, due to this method needs 
several relay feedback tests and only for fixed structure low order models, so it couldn’t meet the actual industrial requirements. 


According to this circumstance, a new technique of modeling identification for process with integral unit model has been raised. 


It can find a main frequency band of the object quickly and analyze frequency response. Then the model with integral unit can 


be identified accurately. The results of simulation show that this method could acquire a great identification effect. 
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般 采 用 继 电 反 馈 法 &53, 该 方法 通过 在 闭环 系统 回路 中 反复 切换 


PID 控制 器 进行 继 电 反 馈 实验 以 采集 多 个 频率 响应 点 进行 计算 。 

在 实际 工业 生产 应 用 中 ， 尤 其 是 在 化 工 工业 以 及 火力 发 电 ”但 是 考虑 到 工业 生产 中 的 经 济 性 、 安 全 性 、 复 杂 性 等 ra, 往往 无 
机 组 等 领域 ， 比 如 水 塔 水 位 控制 、 气 体 存储 中 四 有 时 会 出 现 含 ”法 在 闭环 系统 中 随意 切换 PID 控制 器 进行 数 次 继 电 反馈 实验 ， 
积分 环节 的 过 程 对 象 。 经 典 的 例子 就 是 锅炉 汽 包 外 ， 锅 炉 的 膨 ” ”并 且 该 方法 只 能 辨识 固定 结构 的 对 象 模型 ,所 以 往往 无 法 使 


胀 会 导致 鼓 级 与 锅炉 给 水 量 之 间 的 传递 函数 发 生 改变 ， 产 4 用 继 电 反 馈 法 进行 含 积分 环节 对 象 的 辨识 。 而 文献 [10,11] 所 提 
个 含 积分 环节 的 对 象 模型 。 这 种 对 象 不 同 于 一 般 单调 、 稳 定 的 ”出 的 根据 实际 运行 的 系统 采集 到 的 数据 的 频 域 辨识 方法 ， 也 只 
通用 工业 模型 ， 现 今 学 界 有 很 多 针对 这 种 过 程 对 象 的 开 环 频 域 。 是 针对 一 般 的 过 程 对 象 ， 对 于 这 类 含 积分 环节 的 特殊 对 象 则 不 


FE 


辨识 方法 ， 如 作 图 法 6， 通 过 给 出 系统 开 环 阶 跃 响应 曲线 的 拐 。 ” 太 适 用 。 

点 和 峰值 点 , 即 作出 y = 0.28* y(eo) 处 曲线 的 切线 进行 近似 计算 ; 对 此 为 了 获得 含 积分 环节 对 象 的 准确 模型 ， 本 文 提 出 了 一 
或 者 累加 法 外 , 通过 将 系统 离散 化 后 , 累加 输出 值 进行 计算 。 但” 种 新 的 频 域 辨识 方法 ， 只 需要 采集 闭环 系统 的 实际 运行 数据 进 
是 在 实际 工业 过 程 对 象 的 行 分 析 ， 即 可 找 出 被 辨识 对 象 的 主要 频率 段 ， 然 后 在 该 频率 段 


模型 辨识 中 ， 考 虑 到 系统 的 安全 性 和 经 济 性 ， 无 法 进行 7 了 频率 特性 分 析 ， 就 可 以 迅速 找 出 对 象 的 多 
环 实验 。 而 已 有 的 对 含有 积分 环节 过 程 对 象 的 模型 闭环 辨识 一 ”个 主要 频率 点 ， 从 而 进行 含 积分 环节 对 象 过 程 的 模型 辨识 。 
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1 ”过 程 对 象 的 模型 辨识 


1.1 频率 响应 特性 
通常 在 实际 工业 生产 中 ， 信 号 从 一 个 稳定 状态 变化 到 另 一 
个 稳定 状态 的 情况 经 常 遇 到 ， 两 个 稳定 状态 之 间 的 信号 包含 了 
极为 重要 的 频 域 特性 。 因 为 一 般 的 系统 响应 信号 f(z) 可 以 分 解 
成 瞬 态 部 分 Af (1) 和 稳 态 部 分 上 (D ， 即 
Af(D)= f(D) -£00 (1) 
于 两 个 稳定 状态 之 间 的 变化 的 信号 通常 不 满足 绝对 可 积 
的 条 件 ， 所 以 考虑 对 等 式 两 边 进行 Laplace 变换 ， 得 
F(s)=AF(s)+F(s) 
2 


= = ArDedr+ 一 一 


假设 在 1=7, 时 刻 ， 


本 


时 刻 的 Af(D) 基本 都 等 于 零 ， 则 对 象 的 输入 输出 信号 的 拉 式 变 
换 表 达 式 为 
U(s) = 六 Ax(Dead+zE) G) 
9 8 
cf Ceat (4) 
9 5 
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图 1 含 积分 环节 对 象 模型 的 输入 输出 信和 号 
到 1 所 示 ， 在 闭环 阶 跃 响应 中 ， 对 象 输入 为 阶 跃 信号 ， 
含 动 态 环节 和 稳 态 环节 这 两 部 分 ， 


刀 


五 
双 


Au(t) 等 于 零 ， 而 输出 结果 


如 果 使 得 s = jw ， 则 对 象 的 传递 函数 可 表示 为 : 
ri i Oe 
U (jw) ul(%) 


可 以 将 式 (5) 的 积分 部 分 分 成 个 长 度 为 Ar 的 子 区 间 ， 展 
开 成 积分 累加 的 形式 ， 即 


fed = 2 Ay(De wd (6) 


其 中 N 为 数据 的 采样 个 数 。 
在 式 (6) 的 基础 上 ， 将 s = jw 代入 式 (5)， 进 一 步 扒 
得 到 对 象 各 个 重要 频率 点 的 频率 响应 特性 : 


导 整 理 ， 


Chinaxi 合作 期 干 


过 程 对 象 模型 a 


C(CJw)s 
Ny (Sa 三 下 @ SA ] 
N 大 
es OA 0 
yt De 0 0 
k=1 eo 一 1 1 
Ay. 1 2 = 
OA jw 


u(%) 


频率 w， 都 可 根据 式 (7) 求 出 该 频率 点 时 


这 样 对 于 任意 给 定 
对 象 的 幅 值 和 相位 。 
1.2 过程 对 象 的 主要 频率 点 

含 积 分 环节 模型 的 Nyquist 图 如 图 2 所 示 ， 假 设 过 程 对 象 
的 模型 趋 于 收敛 , 由 于 一 般 系 统 通常 工作 在 中 、 低频 段 外 , 因此 
我 们 只 需要 选取 曲线 的 中 、 低 频段 作为 主要 频率 段 ， 即 主要 频 
率 段 位 于 第 工 II 象限 , 也 就 是 位 于 虚 轴 左 半 平 面 的 频率 点 5 站 。 
假设 频率 为 零 作为 频率 段 的 起 始 频 率 , 临 界 频率 作为 终点 频率 ， 


用 和 迭 代 公式 计算 
O11 = 0, (r+9,) (8) 
9, — Pi 
p=arg[G(j)] 09) 
其 中 :ww 和 gp 均 为 零 ，w 取 一 个 非常 小 的 数 ， 例 如 103 
等 ，g 通过 式 (9) 进 行 计 算得 到 ， 其 中 的 G(jw,) 由 式 (7) 可 得 。 


经 过 数 次 的 反复 迭代 以 后 ， 求 出 对 象 的 临界 频率 w (对 应 相 角 
为 -x )。 此 时 ，(0, wo) 就 是 想 要 求 得 的 主要 频率 段 。 


Nyquist Diagram 


Imaginary Aixs 


Real Axis 


图 2 积分 加 纯 灌 后 模型 的 Nyquist 图 
1.3 传递 函数 模型 

在 根据 2.2 节 方 法 所 求 的 频率 段 (w, wo,) 内 取 M 个 频率 点 
来 求解 最 终 的 过 程 对 象 的 传递 函数 模型 。 在 实际 工业 生产 中 ， 
工业 对 象 的 系统 模型 都 是 高 阶 和 复杂 的 ， 以 至 于 即使 我 们 获得 
其 精确 地 模型 也 很 难 设 计 出 合适 的 控制 器 ， 因 此 我 们 采用 下 式 
(10) 二 阶 带 积 分 模型 对 实际 系统 过 程 进行 逼近 。 

G(s)= i 7 (10) 

式 (10) 的 参数 ab 以 及 工 可 以 通过 将 @ (n=1,2,3.…,M) 代 入 

G (jw) 与 实际 频率 响应 点 G( jw) 进行 拟 合 来 解 得 ， 即 
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b 
G(i 二 -Lon 
U0 ey 
/= M 
其 中 ， 在 拟 合 时 需要 用 到 G( jw) 的 幅 值 条 件 和 相位 条 件 ， 即 
@,|G(CoO + (12) 
w |G(jo| a’ -三 =0 
—arg[G(j%,)]—~— 
2 (13) 
tan-!( 人 ) = 中 ,上 
a 
其 中 =1,2,3,...,M ， 式 (12) 可 以 用 式 (14) 的 和 矩阵 表示 。 
DO=T (14) 
olGGo) olGGo) 
中 二 oo |GCo ooGCoo) (15) 
wlG(jo) wlG(jo)| 
1 1 
1 0 pb? 
二 ，0= 16 
r-l!l,g 四 (10) 
1 | 
式 (14) 中 的 g 可 以 通过 最 小 二 乘 算法 求 出 
0O=(Oro) PT (17) 
式 (16) 可 以 计算 出 传递 函数 模型 的 各 个 参数 
1 
加 m 
中 |, 民 
1Q 
如 果 系 统 的 输出 随 着 输入 信和 号 的 增加 而 增加 , 则 为 正 对 象 ， 
qa,b 取 正 ; 反之 ， 则 取 负 。 
在 模型 参数 a,p 确定 以 后 , 纯 滞 后 工 则 可 以 通过 G(jw) 的 
相位 条 件 和 最 小 二 乘法 求 出 。 
-as[Gual- 二 -am( 包 | 
四 2 a 
， A ff 
有 加 ~arg[G(jo,)] -= -tan (2 (19) 
O, 
—arg[G(jw,)]— tan-! (和 
L 2 a /J 
2 ”仿真 实验 
为 了 检验 本 方法 的 准确 性 和 可 靠 性 ， 本 文 在 仿真 实验 中 引 
入 信 噪 比率 (SNR) [12,13,14,15] 
mean(abs(noise)) (20) 
mean(abs(signal)) 
于 传统 的 继 电 反 馈 法 无 法 辨识 高 阶 对 象 模 型 ， 而 实际 工 
业 中 需要 进行 高 阶 模型 的 辨识 ， 所 以 采 | 高 阶 模型 作为 对 象 模 
设 定 仿真 对 象 一 的 模型 16 为 
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张 ”， 迁 ， 等 : 含 积分 环节 芭 册 恨 呈 衣 红 相 的 频 域 辨识 方法 
1 
= 21 
s(s+1)s 人 
其 开 环 阶 跃 响应 曲线 如 图 1 所 示 。 
_ Nyquist Diagram 
10 + T r 
0 
-10 
-20 
四 30 
三 
-40 
-50 
-60 1 1 1 1 1 1 
-9 -6 -5 二 -3 -2 -1 0 1 
Real Axis 
图 3 对 象 1 在 SNR=20% 环 境 下 辨识 模型 的 Nyquist 图 
图 3 为 在 SNR=20% 环 境 下 的 实际 过 程 对 象 一 和 采用 本 文 
提出 的 方法 所 得 到 模型 的 Nyquist 曲线 。 图 中 红色 实 线 为 实际 


对 象 的 Nyquist 图 ， 蓝 色 实 线 为 本 文 方法 辨识 出 的 模型 的 
Nyquist 图 ， 黑 色 电 话 线 为 一 般 二 阶 辨识 方法 辨识 出 的 模型 的 


Nyquist 图 ， 
最 终 辨识 出 的 主要 频率 范 


x 为 对 象 的 频率 响应 点 。 


围 为 (0, 0.1852) ， 在 此 区 间 选 取 


十 个 频率 点 ， 最 终 辨 识 出 的 对 象 模 型 为 

0.3341 si1w6 
(OY 22 
(3) s(s+0.3336) -9 

设 定 对 象 二 的 模型 为 积分 加 纯 灌 后 模型 

1 
-+ 23 
OD a we 


Nyquist Diagram 


Imaginary Axis 
i 


1 
| 
| 
9 小 ， 
2 


h 1 
-1 -0.8 


图 4 对 象 2 在 SNR=20% 环 境 下 辨识 模型 的 Nyquist 图 
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图 4 为 在 SNR=20% 环 境 下 的 实际 过 
提出 的 方法 所 得 到 模型 的 Nyquist 


十 程 对 象 1 和 采用 本 文 
线 。 图 中 红色 实 线 为 实际 


对 象 的 Nyquist 图 ， 蓝 色 实 线 为 本 文 方法 辨识 出 的 模型 的 


Nyquist 图 ， 
Nyquist 图 ， 


黑色 电话 线 为 一 般 二 阶 辨识 方法 辨识 出 的 模型 的 
x 为 对 象 的 频率 响应 点 。 


最 终 辨 识 出 的 主要 频率 范 
十 个 频率 点 ， 最 终 辨识 


弹 为 (0, 0.5243) ， 在 此 区 间 选 取 


出 的 对 象 模型 为 


录用 稿 
0.6401 ,is 
a 24 
°09 s(s+0.2143) 人 
上 述 两 种 过 程 对 象 的 辩 识 结果 可 以 看 出 ， 本 文 提出 的 频 
或 方法 辨识 出 的 模型 精度 很 高 ， 与 实际 对 象 模型 的 拟 合 程度 也 


很 好 。 采 用 本 文 方法 能 够 根据 采集 到 的 系统 运行 数据 ， 提 炼 出 
所 需 的 重要 频率 段 , 可 以 很 好 的 将 噪声 的 干扰 消除 。 由 此 可 见 

即使 在 噪声 环境 下 ， 本 文 的 辨识 方法 也 具有 很 好 的 本 方法 很 好 
的 鲁 棒 性 和 实用 性 ， 能 够 精确 的 辨识 出 含 积分 环节 的 过 程 对 象 
的 模型 ， 且 本 文 方法 适用 于 广泛 的 含 积分 环节 的 过 程 对 象 。 因 
此 本 文 的 辨识 方法 可 用 于 实际 的 工业 生产 中 。 


;证 


3 ”结束 语 


本 文 提 出 了 一 种 含 积分 过 程 对 象 的 频 域 辨识 方法 ， 不 需要 
任何 的 先 验 条 件 ， 只 需要 根据 系统 正常 运行 产生 的 数据 获得 对 
象 的 主要 频率 段 ， 求 取 重 要 频 域 点 ， 最 终 通过 幅 频 特性 来 确定 
最 终 的 模型 ， 仿 真实 验 表 明 即 使 在 噪声 干扰 下 ， 该 辨识 方法 对 


含 积分 环节 的 模型 仍 具有 良好 的 辨识 效果 ， 弥 补 了 传统 继 电 反 
馈 法 只 能 辨识 固定 的 低 阶 模型 ， 且 需要 大 量 实验 的 不 足 。 并 且 
于 本 文 方法 具有 很 好 的 鲁 棒 性 和 准确 性 ， 适 用 于 广泛 的 含 积 


分 环节 过 程 对 象 ， 所 以 本 文 方法 可 以 用 于 实际 的 工业 对 象 辨识 
中 。 
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